IMBI M ORIGINAL RESEARCH ARTICLE 

inflHfll HVI4"\I f\f*\f published: 21 October 2013 

IIVNVIUIMULUljY doi: 10.3389/fimmu. 2013. 00342 

Association analysis of ILWJNF-a, and IL23R-IL12RB2 
SNPs with Behget's disease risk in Western Algeria 

Ouahiba Khaib Dit Naib\ Mourad Aribi 1 - 2 *, Aicha ldder 3 ,Amel Chiali 4 , Hakim Sairi 4 , Isabelle Touitou 5 - 6 - 7 , 
Gerard Lefranc 8 and Mouna Barat-Houari 5 

' Laboratory of Applied Molecular Biology and Immunology. Department of Biology, Abou-Bekr Belkaid University Tlemcen, Algeria 

2 Departments of Pharmacy and Oral Medicine, Abou-Bekr Belkaid University, Tlemcen, Algeria 

3 Etablissement Hospitaller Specialise en Ophtalmologie, Clinique Hamou Boutlelis, Oran, Algeria 
* Service de Dermatologie, Centre Hospitalo-Universitaire d'Oran, Oran, Algeria 

5 Unite Medicale des Maladies Auto-inflammatoires, Departement de Genetique, CHRU, Montpellier, France 

6 Universite Montpellier 1, Montpellier, France 

1 Genetique des Maladies Auto-inflammatoires etdes Osteoarthropathies chroniques, INSEFIM U844, Montpellier, France 

s Laboratoire d'lmmunogenetique Moleculaire, Institut de Genetique Humaine, CNFtS UPR 1142, et Universite Montpellier 2, Montpellier, France 




Edited by- 
Giuseppe Alvise Ramirez, Universita 
Vita-Salute San Raffaele, Italy 

Reviewed by: 

Fulvio DAcquisto, Queen Mary 
University of London, UK 
Angelo A. Manfred!, Universita 
Vita-Salute San Raffaele, Italy 

"Correspondence: 

Mourad Aribi, Laboratory of Applied 
Molecular Biology and Immunology, 
Abou-Bekr Belkaid University of 
Tlemcen, Rocade # 2, Tlemcen, TLM 
13000, Algeria 
e-mail: m_aribi@mail. 
univ-tlemcen. dz; m_aribi@yahoo. fr 



Objective: We have conducted the first study of the association of interleukin (IU-10, tumor 
necrosis factor alpha (TNF-a), and IL23R-IL 12RB2 region single nucleotide polymorphisms 
(SNPs) with Behget's disease (BD) in Western Algeria. 

Methods: A total of 51 BD patients and 96 unrelated controls from West region of Algeria 
were genotyped by direct sequencing for 11 SNPs including 2 SNPs from the IL10 promoter 
[c.-819T> C (rs1 800871), c.-592A> C (rs1 800872)], 6 SNPs from the TNF-a promoter [c- 
1211T>C (rs1799964), c.-1043C>A (rs1800630), c.-1037C>T (rs1799724), c.-556G>A 
(rs1800750), C.-488G > A (rs1 800629), and C.-418G > A (rs361525)], and 3 SNPs from the 
IL23R-IL 12RB2 region [g.67747415A > C (rs12119179), g.67740092G > A (rs11209032), and 
g.67760140T> C (rs924080)]. 

Results:The minor alleles c.-819Tand c.-592Awere significantly associated with BD [odds 
ratio (OR) = 2.18; 95% confidence interval (CI) 1.28-3.73, p = 0.003]; whereas, there was 
weaker association between TNF-a promoter SNPs or IL23R-IL12RB2 region and disease 
risk. 

Conclusion: Unlike the TNF-a and the IL23R-IL12RB2 region SNPs, the two IL10 SNPs 
were strongly associated with BD.The -81 9T, and -592A alleles and the -819TT, -819CT, and 
-592AA and -592CA genotypes seem to be highly involved in the risk of developing of BD 
in the population of Western Algeria. 

Keywords: Behget's disease, genetic association, IL10, TNF-a, IL23R-IL12RB2, single nucleotide polymorphism. 
Western Algeria 



INTRODUCTION 

Behcet's disease (BD) is a systemic inflammatory multifactorial 
disease (1). It is characterized by recurrent episodes of oral and 
genital ulceration, skin, and ocular lesions (2). The disease is now 
recognized as a systemic vasculitis, given that it can affect other 
tissues and organs including blood vessels, the digestive tract, and 
nervous system (3, 4). 

The etiology of BD is not fully elucidated. There is a hypoth- 
esis that a pathogenic autoimmune process of BD is triggered 
by an infectious or environmental agent, in individuals geneti- 
cally predisposed (5). The most strongly genetic factor associated 
with the disease is HLA B51 (5, 6). This association was initially 
described in 1973 (7) and subsequently confirmed in different 
ethnic groups (6, 8, 9). 

However, the association between BD and HLA B51 represents 
only 20% in the siblings of patients with the disease (10), and 
50% of cases with BD are negative for this allele (8, 9). These 



observations suggest the existence of other risk factors outside 
the HLA region. In fact, several recent genome-wide association 
studies have identified additional and new genomic regions that 
predispose to the disease (11-14). 

The BD causes inflammation and chronic immune activation 
within small blood vessels (15-17). The site of inflammation is 
usually characterized by infiltration of immune cells as well as by 
highly elevated levels of different cytokines (18). IL-10 is one of 
the most important cytokines that has been observed at increased 
level in the serum and active lesions of BD patients (19-21). TNF- 
a and soluble TNF-a receptors are also elevated in the sera of 
patients with BD (22-24). Both IL-10 and TNF-a cytokines have 
been shown to play an important role in the immunopathol- 
ogy of autoimmune diseases (25-28), and an opposite roles in 
the inflammatory responses (29, 30). An autoregulatory loop 
appears to exist in whereby TNF-a induces IL-10 production, 
which ultimately reduces TNF-a synthesis (31, 32). It has been 
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suggested that the increase of IL- 10 may down-regulate the expres- 
sion of NO, prompting the protective role of elevation of IL-10 
(33). Additionally, treatment with anti-TNF-a monoclonal anti- 
bodies has resulted in improvement of various manifestations 
ofBD(34,35). 

IL-10 and TNF-a production may be regulated at the tran- 
scriptional level. Thus, several single nucleotide polymorphisms 
(SNPs) at the promoter of IL10 and TNF-a gene have been shown 
to be associated with changes in the expression levels of IL-10 and 
TNF-a production (36, 37). On the other hand, numerous recent 
studies have demonstrated an association between BD and several 
IL10 (13, 14, 38, 39) and TNF-a (39-43) SNPs in different eth- 
nic groups. However, to date, there are no analogous or identical 
investigations in Algeria. 

Two others cytokines, IL-23 and IL-12, may play an impor- 
tant role in BD pathogenesis; their levels are elevated in BD 
patients (44-46). IL-23 drives and promotes the development 
of a unique T-helper cell population that produces IL-17, 
Thl7 cells. These IL-23-driven Thl7 cells are highly patho- 
genic and elicit IL-17-dependent inflammation in autoimmune 
diseases (47). IL-12, a heterodimeric cytokine, is of crucial 
relevance to cell-mediated immunity and Thl differentiation 
(48). This cytokine exerts its biological effects via binding to 
a heterodimeric receptor consisting of IL12RB2 and IL-12RB1 
subunits. 

The effect of IL-23 and IL-12 is mediated through the IL- 
23 and the IL-12 receptor (IL-23R, IL-12RB1). The genes that 
encode these receptors are adjacent on chromosome lp3 1; a GWAS 
Studies revealed that IL23R-IL12RB2 region is associated with 
BD (13, 14). Nevertheless, its association in the pathogenesis of 
BD remains to be confirmed in different ethnic groups. In this 
context, we examined genetic association for 11 SNPs in IL10, 
TNF-a, and IL23R-IL12RB2 candidate genes with BD in Western 
Algeria. 

MATERIALS AND METHODS 
PATIENTS AND SUBJECTS 

Fifty-one (51) unrelated BD patients and age- and sex-matched 96 
healthy controls originate from the Western Algeria were recruited 
for a case-control study at the Oran Ophthalmic Hamou Boutlelis 
Hospital, the Department of Dermatology of Oran Medical Cen- 
tre University, and the Oran Blood Transfusion Centre (Algeria). 



Among the 51 patients, 11 DNA belonging to Algerian origin, 
were selected from the biobank DNA for Genetics Laboratory 
of Autoinfiammatory Diseases, Arnaud de Villeneuve Hospital, 
Montpellier (France). 

Consent was signed by each participant or participant's parent 
or legal guardian if entrant is a minor, under the Rules of Ethics 
and Professional Conduct. Patient characteristics were recorded 
using a questionnaire. The diagnosis of patients was based espe- 
cially on the criteria proposed in 1990 (49). The control group was 
composed of healthy subjects without a family history of autoin- 
fiammatory diseases, and selected from the same population. This 
work was approved by the Institutional Ethics Board of Tlemcen 
Abou-Bekr Belkaid University. 

GENOTYPING 

Each DNA was genotyped for 11 SNPs, including two 
IL10 promoter SNPs [c.-819T>C (rsl800871), c.-592A>C 
(rsl800872)], six SNPs from the TNF-a promoter [c- 
1211T > C (rsl799964), C.-1043C > A (rsl800630), C.-1037C > T 
(rsl799724), C.-556G > A (rsl800750), C.-488G > A (rsl800629), 
and c.-418G>A (rs361525)], and three SNPs from the IL23R- 
IL12RB2 region [g.67747415A > C (rsl21 19179), g.67740092G > A 
(rsl 1209032), and g.67760140T > C (rs924080) ] . 

Genotyping was performed at the Laboratory of Genetics 
of Autoinfiammatory Diseases, Arnaud de Villeneuve Hospital, 
Montpellier (France). Genomic DNA was isolated from periph- 
eral blood, drawed on EDTA anti-coagulant, using QIAamp DNA 
Blood Kits (Qiagen, Valencia, CA, USA). The DNA samples were 
then dosed by spectrophotometry ND-1000 (Nano Drop Tech- 
nologies, Wilmington, DE, USA) at 260 and 280 nm. The DNA 
concentration and ratio OD260/OD280 were estimated for each 
sample (50). 

The DNA samples were subsequently amplified in a Applied 
Biosystems Thermocycler (Applied Biosystems, Foster City, CA, 
USA) in a 15 \xL reaction volume containing 50 ng DNA, 
2X Promega PCR Master Mix, and 25 u,M of each primer 
(Table 1). The PCR programs were as follows: after a denat- 
uration phase of 15min at 95°C, the samples were subjected 
to 35 amplification cycles followed by a final elongation step 
of 7 min at 72°C. Each cycle comprises 30 s denaturation at 
95°C, 30 s of primer annealing at 60°C, and 1 min extension 
at 72°C. 



Table 1 | Primers sequence and length product. 



Loci 


SNPs 


Forward primer 


Reverse primer 


Product length (bp) 


IL10 


rs1800871 rs1800872 


TTAG ACTC C AG C C AC AG AAG C 


GGGGGACCCAATTATTTCTC 


597 


TNF-a 


rs1799964 rs1800630 rs1799724 


GTGTGTCTGGGAGTGAGAACTTC 


CTTCTTTCATTCTGACCCGG 


570 




rs1800750 rs1800629 rs611525 


CTCAGGACTCAACACAGCTTTTC 


GAAAGAATCATTCAACCAGCG 


438 


IL23R-IL12RB2 region 


rs11209032 


GGAGTTAAACCTCTTGCTATCCTG 


GATGCACAATGAGTTGATAAGG 


164 




rs12119179 


TACCC AG G G CATTCAG CTAC 


G CTTG AGCTCCTG GATC AAG 


701 




rs924080 


GCACGTATGCC I I I I I GCATA 


ATTTGAATGTGCCTTGGCAT 


364 



bp, base pair; IL, interleukin; rs, reference SNP; SNP single nucleotide polymorphism; TNF, tumor necrosis factor. 
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After checking the quality and size of the PCR products by 
agarose gel (1.5%) electrophoresis, SNPs genotyping was per- 
formed by direct sequencing using the BigDye Terminator ver- 
sion 3.1 (BDT v3.1) Cycle Sequencing Kit, followed by capillary 
electrophoresis on an ABI 3100XL Genetic Analyzer, according to 
the manufacturer's recommendations (Applied Biosystems, Foster 
City.CA, USA) (Figure 1). 



STATISTICAL ANALYSIS 

Comparisons of allele and genotype frequencies between groups 
(patients versus control subjects, and between the patient's groups 
according to different clinical features) were performed using the 
Chi-square or Fisher's exact tests. The association analysis was 
carried out by Odds ratio (OR) and corresponding 95% confi- 
dence interval (95% CI). Statistical analyses were performed using 



rs1800871 



lllllll Illllllllllllllllllllllllll 

tcagttggcactggtgtacccttgtacaggtgatgtaa; 



lllllllllllllll ll Illllllllllllllllllllllllll 

STCTCTGTGCCTCAGTTTGCTCACTATAAAATAGAGACGGTAGG 




llllllllllllllllllllllllllllllllllllll 

TCAGTTGGCACTGGTGTACCCTTGTACAGGTGATGTAAfc 



llllllllllllllllllllllllllllllllllllllllllll 

TCTCTGTGCCTCAGTTTGCTCACTATAAAATAGAGACGGTAGG 




llllllllllllllllllllllllllllllllllllll 

TCAGTTGGCACTGGTGTACCCTTGTACAGGTGATGTAAE 



lull I. Iilllll 1 1 1 lllllll iilllll lliilli 

TCAGTTgGCACTGGTGTACCCTTGTACAGGTGATGTAA 



lllllllllllllllllllllllllllllllllllllllllll 

flTCTCTGTGCCTCAGTTTGCTCACTATAAAATAGAGACGGTAGG 




llllllllllllllllllllllllllllllllllllllllllll 

1TCTCTGTGCCTCAGTTTGCTCACTATAAAATAGAGACGGTAGG 




rs1800872 

Mtlllllllllllililllliiii.iillliiiiiilllilillllllllllll.llllll.ini, 

CCTGGAACACATCCTGTGACCCCGCyTGT::TGTAGGAAGCCAGTCTCTGGAAAGwAAAATGGAAGG( 




■III llllllllla III ■■■■ .III 

CCTGGAACACATCCTGTGACCCCGCYTGT 



. ■IIIIliillllllllllll.llli i.l Hi 

= |CTGTAGGAAGCCAGTCTCTGGAAAGwAAAATGgAAGG( 



Illllllllllllllllllllllll.nl 

CCTGGAACACATCCTGTGACCCCGCyTGT 




lllllllllllllllllllllll I ; I I I I I I . I I I I I 

5:TGTAGGAAGCCAGTCTCTGGAAAGwAAAATGGAAGG( 



lllllllllllllllllllllllli.lll 

CCTGGAACACATCCTGTGACCCCGCyTGT 





llillliiillllllllllii.iiliil.ini 

AGGAAGCCAGTCTCTGGAAAGwAAAATGGAAGGI 



FIGURE 1 | Electropherogram of rs1800871 and rs1800872. rs, reference SNP; SNP single nucleotide polymorphism. 
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GraphPad Prism Version 5.04 (GraphPad Software, Inc., La Jolla, 
CA, USA) and Epi Info 2000 Version 1.0 for Windows (Epi Info, 
Atlanta, GA, USA) software. 

RESULTS 

Table 2 shows the description of the clinical characteristics of the 
patients with BD of the current study. The mean age (±SD) of the 
patients at disease onset was 26 ± 1 1 years. Predominant lesions 
were oral ulcers ( 100%), cutaneous lesions (86.27%), genital ulcers 
(82.35%), eye lesions (62.74%), and arthritis (58.82%). 

The distribution of alleles and genotypes frequencies of 
IL10 promoter SNPs c.-819C>T (rsl800871) and c.-592C>A 
(rs 1800872) showed that the two SNPs were in total linkage dise- 
quilibrium in our sample. For this, reason the results of one SNP 
C.-819C > T will be considered (Table 3). 

The allele frequencies were significantly different in patients 
compared to controls. As indicated in Table 3, the frequencies of 
C.-819T (rsl800871) allele, and of the -819TT,-819CT (rsl800871) 
genotypes were significantly increased in patients than in controls 
(p = 0.003 and p = 0.005, respectively). Additionally, these SNPs 



Table 2 | Clinical and demographic features of the Behcet patients of 
the current study. 



Characteristics 

Mean age at disease onset ±SD (year) 

Sex ratio M/F (%, n) 

Oral ulcers (%, n) 

Genital ulcers (%, n) 

Cutaneous lesions (%, n) 

Eye lesions (%, n) 

Neurological symptoms (%, n) 

Venous thrombosis (%, n) 

Arthritis (%, n) 

Multiplex family (%, n) 

pediatric case (%, n) 

consanguinity (%, n) 



Frequency [n = 51) 



26±11 
56.9/43.1 
100 (51) 

82.4 (42) 
86.3 (44) 

62.7 (32) 
35.3 (18) 

25.5 (13) 

58.8 (30) 
35.3 (18) 

19.6 (10) 
43.1 (22) 



(29/22) 



was significantly associated with the disease (c.-819T; OR = 2.18, 
95% CI 1.28-3.73, p < 0.01; -819TT and -819CT, OR = 2.17, 95% 
CI 1.01-4.69, p < 0.05) (Figure 2). 

A subset analysis was performed to examine the difference in 
allele frequencies in clinical subsets of BD (Table 4). We observed 
a significant association between C.-819T and all classes; never- 
theless, the association was slightly lower for the ocular lesion 
(OR= 1.55, 95% CI 0.81-2.96, p > 0.05). Additionally, the associ- 
ation was more significant for the Genital ulcers (OR = 2.21; 95% 
CI 1.29-4.04, £ = 0.002). 

We reported in Tables 5 and 6 that all IL23R-IL12RB2 SNPs 
alleles and genotypes, respectively, were not significantly associated 
with the disease (OR > 1, p > 0.05). The minor allele frequencies 
were different in the two groups, but this difference did not reach 
statistical significance (p > 0.05). 

As indicated in Tables 5 and 7, alleles and genotypes of the 
TNF-a polymorphisms display similar distributions in patients 
and controls (p > 0.05). Except for c- 1037T and C.-488A all others 
TNF-a alleles were not associated with BD (OR < 1). 



SD, standard deviation. 



2.182 
2.18 

.2 2178 
2 2.176 
in 

2 2.174 
° 2.172 
2.17 
2.168 



C.-819T 



-819TT 
-819CT 
— • 



0"h 



0 1 2 3 4 5 
95% CI 

FIGURE 2 | Odds ratios for associations between IL10 C.-819C >T, and 
-819TT, -819CT (rs1800871) with Behcet's disease. The two SNPs 
C.-819C >T and C.-592C > A are in total linkage disequilibrium in our sample; 
therefore, the results of only one SNP C.-819C >T is considered. CI, 
confidence interval; rs, reference SNP; SNP single nucleotide 
polymorphism. 



Table 3 | Allelic and genotypic frequencies of rs1800871 variant in BD Table 4 | Association analysis of clinical subclasses with IL10 C.-819T 
patients and controls. SNP in patients with Behcet's disease. 



Alleles and genotypes 


Frequencies (%) 


P 


Clinical subset 


OR 




95% CI 


P 


Controls (n = 96) 


Cases (n = 51) 


LL 


UL 


C 


141 (73.4) 


57 (55.9) 


0.003** 


Eye disease 


1.55 


0.81 


2.96 


0.152 


T 


51 (26.6) 


45 (44.1) 




Genital ulcers 


2.21 


1.29 


4.04 


0.002** 


CC 


50 (52.1) 


17 (33.3) 


0.005** 


Skin lesions 


2.07 


1.17 


3.68 


0.007** 


CT 


41 (42.7) 


23 (45.1) 




Arthritis-arthralgia 


2.06 


1.11 


3.81 


0.013* 


TT 


5(5.2) 


11 (21.6) 




Neurologic signs 


2.5 


1.18 


5.32 


0.009** 



BD, Behget's disease; SNR single nucleotide polymorphism; rs, reference SNR CI, confidence interval; LL, lower limit; UR upper limit; OR, odds ratio; SNR single 
**p<0.01. nucleotide polymorphism. *p<0.05, **p<0.01. 
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Table 5 | IL23R-IL12RB2 and TNF-a allelic frequencies in patients with Behcet's disease. 

Loci SNPs Alleles frequency OR (95% CI) p 

(proportion, %) 



IL23R-IL12RB2 rs12119179 (g. 6774741 5A > C) A C 1.24(0.72-2.13) 0.415 

Patients 67 (65.7) 35 (34.3) 

Controls 135 (70.3) 57 (29.7) 

rs1 1 209032 (g.67740092G> A) G A 1.18(0.68-2.03) 0.530 

Patients 67 (65.7) 35 (34.3) 

Controls 133 (69.3) 59 (30.7) 

rs924080(g. 677601 40T>C) T C 1.45(0.87-2.41) 0.133 

Patients 54(52.9) 48(47.1) 

Controls 84 (43.8) 108 (56.3) 

TNF-a rs1799964(c.-1211T>C) T C 0.92(0.52-1.61) 0.751 

Patients 73 (71.6) 29 (28.4) 

Controls 134 (69.8) 58 (30.2) 

rs1800630 (c.-1043C>A) C A 0.83(0.41-1.67) 0.584 

Patients 86(84.3) 16(15.7) 

Controls 157(81.8) 35(18.2) 

rs17999724 (c.-1037C>T) C T 1.01(0.35-2.84) 0.976 

Patients 95(93.1) 7(6.9) 

Controls 179(93.3) 13(6.7) 

rs1800750 (c.-556G>A) G A 0.61(0.14-2.09) 0.402 

Patients 98(96.1) 4(3.9) 

Controls 180(93.7) 12(6.3) 

rs1 800629 (c.-488G> A) G A 1.12(0.56-2.26) 0.726 

Patients 85(83.3) 17(16.7) 

Controls 163(84.9) 29(15.1) 

rs361525(c.-418G>A) G A 0.93(0.39-2.21) 0.869 

Patients 92 (90.2) 10 (9.8) 

Controls 172 (89.6) 20 (10.4) 



CI, confidence interval; IL23R-IL 12RB2, region of the interleukin-23 receptor and interleukin-12 receptor beta 2 chain-encoding gene; OR, odds ratio; rs, reference 
SNP; SNR single nucleotide polymorphism; TNF, tumor necrosis factor. 



Table 6 |The distribution of IL23R-IL12RB2 genotypes in patients with Behcet disease. 

SNPs Genotype distribution p-Value OR (95% CI) p-value MAF (%) OR MAF (95% CI) p-value 

(frequency, %) 



rs12119179 (g.67747415A > C) 


AA 


AC 


CC 


0.522 


1.13 (0.54-2.35) 0.734 


C 


1.24 (0.58-2.67) 0.580 


Patients 


24 (47.1) 


19 (37.2) 


8 (15.7) 






34.31 




Controls 


48 (50) 


39 (40.6) 


9 (9.4) 






29.69 




rs1 1209032 (g.67740092G > A) 


GG 


GA 


AA 


0.721 


1.12 (0.57-2.21) 0.744 


A 


1.18 (0.71-1.97) 0.560 


Patients 


23 (45.1) 


21 (41.2) 


7 (13.7) 






34.31 




Controls 


46 (47.9) 


41 (42.7) 


9 (9.4) 






30.73 




rs924080 (g. 677601 40T> C) 


CC 


CT 


TT 


0.327 


1.46 (0.64-3.35) 0.326 


C 


1.58 (0.71-3.54) 0.264 


Patients 


13 (25.5) 


22 (43.1) 


16 (31.4) 






47.1 




Controls 


32 (33.4) 


44 (45.8) 


20 (20.8) 






56.25 





ORs were calculated for the minor versus major alleles. CI, confidence interval; IL23R-IL12RB2, interleukin-23 receptor and interleukin-12 receptor beta 2 
chain-encoding gene; MAR minor allele frequency; OR MAf , OR of the MAF; OR, odds ratio; rs, reference SNP; SNR single nucleotide polymorphism. 
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Table 7 |The distribution of TNF-a genotypes in patients with Behcet disease. 

SNPs Genotype distribution p-Value OR (95% CI) p-value MAF (%) OR MAF (95% CI) p-value 

(frequency, %) 



rs 1 799964 (c.-1 2 1 1 T > C) 


TT 


TC 


CC 


0.114 


0.67 (0.32-1.4) 0.249 


c 


0.8 (0.37-1.75) 0.580 


Patients 


29 (56.9) 


15 (29.4) 


7 (13.7) 






28.43 




Controls 


45 (46.9) 


44 (45.8) 


7 (7.29) 






30.21 




rs1800630 (c.-1043C>A) 


CC 


CA 


AA 


0.826 


0.79 (0.35-1.78) 0.544 


A 


0.78 (0.31-1.97) 0.600 


Patients 


37 (72.5) 


12 (23.5) 


2 (4) 






15.69 




Controls 


65 (67.7) 


27 (28.1) 


4 (4.2) 






18.23 




rs17999724 (C.-1037C >T) 


CC 


CT 


TT 


ND 


1 23 (0 4-3 74) 0 690 


T 


1 26 (0 34-4 69) 0 733 


Patients 


44 (86.3) 


7 (13.7) 


0 (0) 






6.86 




Controls 


83 (86.5) 


13 (13.5) 


0 (0) 






6.77 




rs1800750 (c.-556G>A) 


GG 


GA 


AA 


ND 


0.6 (0.13-2.12) 0.388 


A 


0.6 (0.12-3.07) 0.536 


Patients 


47 (92.2) 


4 (7.8) 


0(0) 






3.92 




Controls 


84 (87.5) 


12 (12.5) 


0(0) 






6.25 




rs1800629 (c.-488G>A) 


GG 


GA 


AA 


0.909 


1.17 (0.52-2.61) 0.680 


A 


1.18 (0.47-2.97) 0.721 


Patients 


35 (68.6) 


15 (29.4) 


1 (2) 






16.67 




Controls 


69 (71.9) 


25 (26) 


2 (2.1) 






15.1 




rs361525 (c.-418G>A) 


GG 


GA 


AA 


0.299 


0.75 (0.28-2.02) 0.541 


A 


0.9 (0.29-2.79) 0.849 


Patients 


43 (84.3) 


6 (11.8) 


2 (3.9) 






9.8 




Controls 


77 (80.2) 


18 (18.8) 


1 (1) 






10.42 





ORs were calculated for the minor versus major alleles. CI, confidence interval; MAF minor allele frequency; ND, not defined; OH, odds ratio; OR MAF , OR of the MAF; 
rs, reference SNP; SNP single nucleotide polymorphism; TNF tumor necrosis factor. 



DISCUSSION 

To date, the etiopathogenesis of BD is not fully elucidated. 
Researches in recent decades have shown the complex role of 
genetic factors in the development of the disease. We analyze the 
association between BD and 11 SNPs in IL10, TNF-a, and IL23R- 
IL12RB2 candidate genes in the Western Algeria population. 

This is the first report demonstrating that the C.-819T and 
C.-592A alleles were associated with BD in Algeria. 

Previous genetic studies have shown a strong association of 
many IL10 variants with BD in different ethnic groups. Recent 
genome-wide association study conducted by Mizuki et al. (14) 
in a Japanese cohort, including 612 individuals with BD and 740 
unaffected individuals controls from different ethnic groups, has 
shown a significant difference between the two groups for five 
IL10 SNPs. The two SNPs analyzed in our study showed an asso- 
ciation (OR = 2.18, 95% CI 1.28-3.73, p<0.01). On the other 
hand, a strong association has also been highlighted for IL10 
rsl5181 11 (OR= 1.45, 95% CI 1.34-1.58) in a genome-wide study 
performed in 311,459 SNPs in 1215 individuals with BD and 
1278 healthy controls from Turkey (13). Our results are consis- 
tent with those obtained by Wallace et al. (38) who studied the 
segregation of two IL10 promoter SNPs, rsl800871 and rsl800896 
in 178 cases and 295 controls from two populations, including 
Arab Middle East and United Kingdom. A strong association of 
the 819T allele has been observed in United Kingdom patients 
(OR 1.5, 95% CI 1.1-2). 

Our results are, however, in disagreement with those of the 
study of Ates et al. (39) in which no significant association was 



revealed by exploring three IL10 SNPs (-1082G > A,rsl800896, c- 
819C > T, rsl800871, and C.-592C > A, rsl800872) in 102 patients 
with BD and 102 controls from Turkey. These conflicting results 
may be explained by ethnic differences. 

Several studies have shown the association of the two (c- 
819 C > T, C.-592C > A) studied SNPs with various inflammatory 
diseases, but also with cancer (51-53), periodontitis (54), and 
docetaxel-induced liver injury (55). These associations suggested 
that the two SNPs can play an important role in the expression 
of IL10. In fact, it has been previously reported that the SNP c- 
819C>T and/or c.-592C>A alleles affect the transcription of 
IL10 (56). Other studies performed on three SNPs at position - 
1082 A > G, -819C > T, and -592C > A in the promoter region of 
the IL10 gene have shown that the expression levels of IL-10 was 
significantly different according to the some haplotypes (57, 58). 
Finally, it has been reported that the disease-associated rsl518111 
allele seems to be associated with low IL-10 mRNA expression 
and protein production (13). Indeed at the SNP rs 1518111 locus, 
the rate of transcript of the G allele is higher than that of the 
A allele and patients homozygous for the A allele produce less 
IL-10 than those who are heterozygous or homozygous for the 
G allele. 

The result of subset analysis suggests that the risk allele might 
predispose to genital ulcers, skin lesions, neurologic signs, and 
arthritis-arthralgia, but weakly to eye complications; the asso- 
ciation was more significant for the genital ulcers (OR = 2.21, 
p = 0.002). Our results are in agreement with those of pre- 
vious study (21) that showed a variable increase in mRNA 
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expression within all BD lesions, including oral and genital ulcers, 
pseudofolliculitis lesions, and lesions at the site of pathergy 
testing. 

Recent GWAS study from Turkey and Japan revealed 
IL12R-IL23RB2 SNPs in association with BD. Three SNPs 
were strongly associated with the disease, including rs924080 
(OR = 1.28, p = 6.69 x 1(T 9 ) (13), rsl2119179 (p = 2.7x 1(T 8 ), 
and rsl495965 (OR= 1.35, p = 1.9 x KT 11 ) (14), but no signifi- 
cant association was found in a Korean cohort. 

Our results showed no significant association between 
BD and rsl2119179, g.67740092G > A (rsll209032), and 
g.67760140T > C (rs924080) SNPs in the IL23R-IL12RB2 region. 
In Iranian study (59), six SNPs in IL23R-IL12RB2 were found 
to be associated with BD; the most significant of which were 
rsl7375018 (OR= 1.51, p= 1.93 x 1(T 6 ), rs7517847 (OR= 1.48, 
p= 1.23 x 1(T 6 ), and rs924080 (OR = 1.29, p = 1.78 x 1(T 5 ). 
Others studies have also identified a strong relationship between 
polymorphisms of IL23R and BD (60-62). These associations may 
suggest an important role of Thl7 cells that express the IL-23R on 
their surface. Kim et al. (62) studied the interaction of specific 
IL17A, IL23R, and STAT4 (signal transducers and activators of 
transcription 4) SNPs in intestinal BD Korean patients; they sug- 
gest that the IL-23/IL-17 axis plays a significant role in disease 
pathogenesis. 

IL-12 has been implicated in the pathogenesis of a multi- 
tude of diverse autoimmune diseases (63, 64). IL12RB2 con- 
stitute a risk factor for primary biliary cirrhosis, with the 
reported top associated SNPs mainly located in intronic sequences 
(65-67). 

The genetic architecture and modularity of human autoim- 
mune diseases is very complex. The functional implications 
of most of these associations are not yet clarified. Iden- 
tify candidate causal SNPs and pathways (ICSN Pathway) 
analysis may act as a powerful guide to further research 
into the functional and immunological ramifications of these 
associations. 

No significant associations were found between BD and stud- 
ied TNF-a polymorphisms. These SNPs have been studied in 
various ethnic groups for possible association with BD. How- 
ever, the allelic and genotypic associations of these studies 
have been contradictory. In Korean patients, TNF-a C.-1043A 
(rsl800630) allele was associated with an increased risk of BD 
(OR=1.4, p = 0.030) (68). However, no significant association 
was found in meta-analysis studies for this SNP (42). Addi- 
tionally, it has been reported a significant associations between 
C.-1037T allele (rsl799724) (OR = 0.76, 95% CI 0.58-0.98), c- 
488G allele (rsl800629) (OR= 1.8, p = 0.010) (68), and C.-418A 
allele (rs361525) (OR= 1.51, 95% CI 1.12-2.04) (42), and BD. 
Moreover, no significant associations were identified with other 
TNF-a promoter polymorphisms, such as C.-1037T (rsl799724), 
c.-488A(rsl800629),and C.-556G > A(rsl800750) alleles withBD 
in Moroccan patients (43 ) . The TNF-a c. - 1 2 1 1 C allele (rs 1 799964) 
presented a significant association with BD in several popu- 
lations, including Turkish (p = 0.023) (69), Korean (p = 0.030, 
OR= 1.4) (68), and UK white Caucasoid population (RR = 2.3, 
p = 0.00004) (40). The frequency of the TNF - a c. - 1 2 1 1 C allele was 



significantly higher in Behcet's patients than in healthy controls 
in Moroccan and Tunisian populations (OR=1.65, p = 0.015; 
OR= 1.68, p = 0.02, respectively) (41, 43) and in meta-analysis 
(OR= 1.35, 95% CI 1.09-1.68) (42). This polymorphism has 
been associated with several extra-intestinal manifestations of 
Crohn's disease, including uveitis, erythema nodosum, and large 
joint arthropathy (70), all of which are known to be asso- 
ciated with BD. Further investigation is necessary to deter- 
mine the functional significance of TNF-a C.-1691042C and 
how it participates in the inflammatory dysregulation associated 
with BD. 

Thus, polymorphisms at positions c- 121 IT > C, C.-1043C > A, 
c.-1037C>T, and c.-488G>A have been associated with 
increased transcriptional activity and production of TNF-a in 
some studies (37, 71), in contrast to others (72-75). The over 
production of TNF-a during the course of BD may result in 
other TNF-a polymorphisms or post-transcriptional mecha- 
nisms. Furthermore, TNF-a production is not only under the 
control of the promoter region of TNF-a, and it may also 
result from complex cis and trans interactions among other 
cytokines. 

TNF-a is encoded in the HLA complex on chromosome 6, a 
region that has long been known to be associated with BD. This 
gene-dense region, presentes a strong linkage disequilibrium (76). 
The association between BD and TNF-a could therefore be a result 
of linkage disequilibrium with alleles within this group. So it will 
be interesting to investigate other genes polymorphism among this 
region in our population. 

In conclusion, we replicate the associations between BD and 
the SNPs from the IL23R-IL12RB2 region and c.-1037C>T 
and c.-488G>A TNF-a promoter SNPs. IL10 promoter SNPs 
(rsl800871 and rsl800872) is strongly associated with BD in the 
population of the Western Algeria. It would be interesting to study 
other SNPs to identify additional associations with BD in the 
studied population. 
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